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Abstract. Let Q be an affine semigroup generating Z'', and fix a finitely generated 
Z'^-graded module M over the semigroup algebra k[Q] for a field k. We provide an al- 
gorithm to compute a minimal Z'^-graded injective resolution of M up to any desired 
cohomological degree. As an application, we derive an algorithm computing the local 
cohomology modules H\{M) supported on any monomial (that is, Z'^-graded) ideal /. 
Since these local cohomology modules are neither finitely generated nor finitely co- 
generated, part of this task is defining a finite data structure to encode them. 
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1. Introduction 

Injective resolutions are fundamental homological objects in commutative algebra. 
For general noetherian rings with arbitrary gradings, however, injective modules are so 
big, and injective resolutions so intractable, that effective computations are never made 
using them. But when the ring in question is an affine semigroup ring of dimension d, 
the natural grading by l/" is substantially better behaved: Z'^-graded injective modules 
can be expressed polyhedrally and are therefore quite explicit. In this paper we provide 
algorithms to compute Z'^-graded injective resolutions over affine semigroups rings. 
Part of this task is finding a finite data structure to express the output. 

As an application, we provide an algorithm to compute the local cohomology, sup- 
ported on an arbitrary monomial ideal, of a finitely generated Z'^-graded module over 
a normal affine semigroup ring. As far as we are aware, this is the first algorithm to 
compute local cohomology for any general class of modules over any class of nonregular 
rings. Our motivation was to make a systematic study of conditions on the support 
ideal and the ambient ring that cause local cohomology to have infinite Bass numbers. 

Date: 15 September 2003. 
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That such infinite behavior occurs only over nonregular rings necessitated our working 
over affine semigroup rings, which seem to be the simplest available singular rings. 

To make our context precise, let Q C be an affine semigroup, that is, a finitely 
generated submonoid of Z'^. We assume that Q is sharp, meaning that Q has no 
units, and that Q generates if as a group. Consider the semigroup algebra k;[(5] = 
©aGQ {-^"} '-'^^^ ^ k. The modules that concern us comprise the category M. of 
Z'^-graded modules H = ©^ggd Ha for which there exists a bound independent of a on 
the dimensions of the graded pieces Ha as vector spaces over k. The injective objects 
in Ai are described in Section |2l and every finitely generated Z'^-graded module lies 
in Ai. Our main theorem concerning injectives is the following. 

Theorem 1.1. Fix a finitely generated l/--graded module M over an affine semigroup 
ring k[Q] and an integer n > 0. The first n stages in a minimal l/--graded injective 
resolution of M can he expressed in a finite, algorithmically computable data structure. 

A more precise version, along with a pointer to the algorithms that do the job, is 
stated in Theorem 14.71 The data structure consists of a list of monomial matrices, 
as we define in Section |21 generalizing those for Q = N'' in jMilOOj . The idea of the 
algorithm in Theorem ll.ll is to do all computations using irreducible resolutions |Mil02j 
as faithful approximations to injective resolutions. Background on irreducible hulls is 
presented in Section |21 the algorithms for working with them constitute Section |21 
The derivation of an algorithm for injective resolutions is then completed in Section^ 

Even more seriously than is the case with injective resolutions, a substantial part of 
building an algorithm to compute local cohomology is finding a finite data structure to 
express the output. Indeed, unlike injectives in our category and in stark contrast 
with the regular case (even without a grading |HS93| |Lyu93| ), the local cohomology 
Hj{M) often has neither a finite generating set nor a finite cogenerating set |Har70t 
IHM03j . This remains true even when M is finitely generated and / C k[(5] is a Z*^- 
graded ideal — that is, generated by monomials. Our solution is to decompose Z'^ into 
tractable regions on which the local cohomology is constant. 

Definition 1.2. Suppose if is a Z'^-graded module over an affine semigroup ring k[Q]. 
A sector partition of H is 

1. a finite partition Z'^ = IJse5 lattice Z"^ into sectors, each of which is 
required to consist of the lattice points in a finite disjoint union of rational polyhe- 
dra defined as intersections of half-spaces for hyperplanes parallel to facets of Q; 

2. a finite-dimensional vector space Hs for each sector S E S, along with isomor- 
phisms Ha — > Hs for all Z'^-graded degrees a E S; and 

3. vector space homomorphisms Hs — > Ht whenever there exist a G S and /3 G T 
satisfying j3 — a G Q, such that for all choices of a and (3, the diagram commutes: 

Ha > Hi3 

i i 

Hs Hx 

Write iS h if to indicate the above sector partition. (The commutativity of the above 
diagram implies immediately that x'^"'^ is the identity, and that x^~-^X"^"'^ = x^~'^.) 
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The finite data structure of a sector partition S \- H, including tlie spaces Hs and 
tlie maps x^~'^, clearly suffice to reconstruct H up to isomorphism. The second half 
of this paper is devoted to computing sector partitions for H when H = H]{M) is a 
local cohomology module. 

Theorem 1.3. For any finitely generated if -graded module M over a normal semi- 
group ring k[(5] and any monomial ideal I , each local cohomology module H]{M) has 
an algorithmically computable sector partition S h H]{M). 

Section El demonstrates how sector partitions arise for the cohomology of any com- 
plex of injectives over a normal semigroup ring. Algorithms for producing these sector 
partitions, particlularly the expressions of sectors as unions of polyhedral sets of lattice 
points, occupy SectionlHl The proof of Theorem ll.3l by expressing local cohomology as 
the cohomology of a complex of injectives (algorithmically computed by Theorem I l.lj) 
in the usual way, occurs in Section [7| That section also treats complexity issues. The 
main thrust is that for fixed dimension d, the running times of our algorithms are all 
polynomial in the Bass numbers of the finitely generated input module M and the 
number of facets of Q, times the usual factor arising from the complexity of Grobner 
basis computation, where it occurs. If d is allowed to vary, then the numbers of 
polyhedra comprising sectors increase exponentially with d. 

Theorem II . 31 allows the computation of many features of local cohomology modules. 
For example, Hilbert series simply record the vector space dimensions in each of the 
finitely many sectors S E S. Our algorithms can actually calculate these dimensions 
without computing the maps in part 3 of Definition II. 21 making it easier to determine 
when (for example) H}{M) is nonzero. Future algorithmic methods (currently open 
problems) include the calculation of associated primes and locations of socle degrees 
(even if there are infinitely many) using a sector partition as input. In particular, 
because of the finiteness of the number of polyhedra partitioning sectors, we believe 
that the socle degrees should lie along polyhedrally describable subsets of Z'^. 

Historical context. There have been a number of recent algorithmic computations 
in local cohomology, such as those by Walther |Wal99j (based on abstract methods 



of Lyub eznik |Lyu93| ), Eisenbu d-Mustata-Stillman |EMS()()j . Miller |Miinnj . Mustafa 
|MusOOj , and Yanagawa |Yan02j . These and related papers fall naturally into a number 
of categories. For instance, the last three deal with Z'^-graded modules over polynomial 
rings in d variables; in particular, they compute local cohomology with support on 
monomial ideals. In contrast, the paper pMSOO works with coarser gradings — but 
still with monomial support, while |Wal99j requires no grading at all. As the gradings 
used become coarser, the papers increasingly depend on Grobner bases: the monomial 
ideal papers require very little (if any) Grobner basis computation; the coarser gradings 
depend heavily on commutative Grobner bases; and the nongraded methods rely on 
noncommutative Grobner bases over the Weyl algebra. 

Regardless of the methods, all of the above papers share one fundamental aspect: 
the base ring is regular (usually a polynomial ring, in the algorithmic setting). The 
reason for restricting to these rings is that local cohomology over them behaves in 
many respects like a finitely generated module, even though it usually fails to be 
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finitely generated. For example, Lyubeznik |Lyu93| and Walther |Wal99j take advan- 
tage of the fact that local cohomology modules over regular rings are finitely generated 
(indeed, holonomic) over the corresponding algebra of differential operators, and that 
the algebra of differential operators of a regular ring is easily presented, at least in 
characteristic zero. 

Generally speaking, our methods lie somewhere between the monomial and coarsely 
graded methods described above, relying on a mix of Grobner bases and integer pro- 
gramming. The principle underlying our computation of injective resolutions is that 
one should attempt to recover entire Z'^-graded modules from their Q-graded parts. 
This idea originated for polynomial rings in |Mil98t IMusOO| IMil OO , was transfered in a 



restricted form to semigroup rings in [YanOlj , and developed generally for semigroup- 



graded noetherian rings in |HM03j . In the present context, the recovery of a module 
from its Q-graded part suggested that we compute injective resolutions via the irre- 
ducible resolutions of |Mil02j . 

Origins of the notion of sector partition can be seen in the Hilbert series formula for 
the local cohomology of canonical modules of normal semigroup rings |Ter99| IYan02j , 
where the cellular homology was constant on large polyhedral regions of li^ . The 
accompanying notion of straight module [yanOll IHM03j abstracted this constancy; 
in fact, our Theorem 15.21 is really a theorem about straight modules as in |HM03t 
Definition 5.1]. In any case, once the injective resolution has been computed using 
irreducible resolutions, the sector partition for local cohomology requires the entire 
Z'^-graded structure of the injective resolution, and not just its Q-graded part. 

Conventions and notation. In addition to the notation introduced thus far, we 
close this Introduction with a note on conventions. The semigroup Q is required to be 
saturated in Sections EHZI because we do not know how to compute sector partitions in 
the unsaturated context f Remark 16. 7|) . Other than the temporary saturation require- 
ment in Section 13.21 the semigroup can be unsaturated in Sections |2}H (Reminders 
of these conventions appear in each section). 

The symbol x" G k[Z'^] denotes a Laurent monomial in the localization k[Z^] of the 
semigroup ring k[Q]. The k- vector space spanned by {x" \ a E T} for a subset T C Z'^ 
will be denoted by k{T}. The k-subalgebra of k[Z'^] will be denoted by k[T]. 

The faces of Q are those subsets minimizing linear functionals on Q. The edges and 
facets are the faces of dimension 1 and codimension 1. To every face F corresponds 
a prime ideal Pp and a quotient affine semigroup ring k[F] = k{F}. 

All modules in this paper are Z'^-graded unless otherwise stated. In particular, 
injective modules (defined in Section^)) are Z'^-graded injective, which means that they 
are usually not injective in the category of all k[(5]-modules. Two subsets S", T C Z'^ 
have the difference set T — S = {P — a \ a & S and /? G T} C Z*^. This allows us to 
write the localization of M along a face F as the module M[ZF] := M ®k[Q] k[Q — F]. 
Homomorphisms N ^ N' of modules are assumed to have Z'^-graded degree 0, so that 
No, N'o, for all a e Z"^. 

We assume in this paper that standard algorithmic calculations with finitely gener- 
ated modules over k[Q] are available. In particular, we assume that the homology of 
any three-term (nonexact) sequence of finitely generated modules can be calculated. 
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as can the submodule annihilated by a prime ideal of k.[Q]. The Z'^-grading only makes 
these computations easier, and the results of all such algorithms are still Z'^-graded. 



In this section the affine semigroup Q need not be saturated. In the Z -graded 
category Ai from the Introduction, the injective modules have simple descriptions. 

Definition 2.1. Let T C Z'^ be closed under addition of elements of —Q, by which we 
mean T — Q <Z T. Then k{T} can be given the structure of a k[(5]-module by setting 



An indecomposable injective is any module of the form k{a + F — Q}, for some 
face F and a G Z'^. 

All such objects are injective in Ai, and every injective object of A4 is isomorphic 
to a finite direct sum of indecomposable injectives |MS03| Chapter 11]. We shall work 
exclusively with objects in Ai. Thus the term "injective module" in the rest of this 
paper will refer to modules of the above type. 

Injectives are infinitely generated. For computations, we therefore work with certain 
finitely generated approximations. A module is called Q-graded if A^ equals its 
Q-graded part Nq := ^^eQ ^ submodule A^ of a module A^' is an essential 
submodule if A^ intersects every nonzero submodule of A^' nontrivially; the inclusion 
A^ "-^ A^' is also called an essential extension. In particular, A^ must be nonzero. 

Definition 2.2. An irreducible sum is a module that can be expressed as the Q- 
graded part Jq of some injective module J. An irreducible hull of a Q-graded 
module A^ is an irreducible sum W along with an essential extension A^ ^ W. 

The existence of unique minimal injective resolutions [MSOS, Corollary 11.35] in- 
cludes the fact that every finitely generated module has an injective hull (that is, an 
inclusion into an injective that is an essential extension) that is unique up to isomor- 
phism. Taking Q-graded parts yields immediately the following lemma. 

Lemma 2.3. Every Q-graded module has an irreducible hull. It is unique up to iso- 
morphism, and isomorphic to the Q-graded part of an injective hull of M. 

We call the modules of Definition 12.21 irreducible sums because of the next lemma, 
which is [Mil02, Lemma 2.2]. An ideal W is called irreducible if W can not be 
expressed as an intersection of two ideals properly containing it. 

Lemma 2.4. A monomial ideal W is irreducible if and only if the Q-graded part of 
some indecomposable injective module J satisfies Jq = W. 

Modules M are usually stored as data structures keeping track of their generators 
and relations — that is, as quotients of free modules. In the context of injective res- 
olutions and local cohomology, storing M as a submodule of an irreducible sum is 
also useful. Our next definition specifies a data structure that precisely describes an 
irreducible sum W. 



2. Effective irreducible hulls 
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Definition 2.5. Effective data for an irreducible sum W = ^{aj + Fj — Q}q 

consist of: 

1. an ordered r-tuple Fi, . . . ,Fr of faces of Q; and 

2. an ordered r-tuple where aj G Z'^/ZFj satisfies Q n {aj + ZFj) ^ (Is. 
An effective vector of degree a G Q is an r-tuple (Ai, . . . , A^) G such that Xj = 
whenever a ^ aj + Fj — Q. Concatenation of the respective face and degree data from 
two effective data yields their direct sum. 

Note that the faces Fj need not be distinct. The condition a G Z'^/ZF takes 
care of the fact that two degrees a and a' off by an element of ZF give the same 
module k{a + F — Q} = k{a' + F — Q}. Usually the a's are recorded as elements 
of Z*^, since the quotient mod ZF can be deduced from the face data. The condition 
Q D {a + ZF) 7^ ensures that k{a + F — Q} has nonzero Q-graded part. The condi- 
tion on the A's simply requires each nonzero component to lie in a nonzero degree of 
the corresponding irreducible summand. 

Definition 2.6. An effective irreducible hull of a Q-graded module M consists of 
effective data for W plus a list of finitely many effective vectors in W generating a 
submodule isomorphic to M. 

An irreducible hull M W is not quite dual to an expression .7-" — > M as a quotient 
of a free module. The generators of JF have as their dual notion the face data Fi, . . . , Fr, 
which as abstract objects associated to M are known as cogenerators. Just as the 
degrees of the generators of J-' need to specified, so must the degree data for the 
cogenerators. However, the notion of effective vector for M as a submodule of W is dual 
not to the notion of relation for M inside J-', but rather to the notion of cogenerator 
for M. Relations for M are, in actuality, dual to the notion of cogenerators for the 
cokernel of M W, which correspond to indecomposable summands in cohomological 
degree 1 of the minimal injective resolution of M; we dub these the correlations of M. 
Thus a presentation of M by generators and relations is dual to a presentation of M 
by cogenerators and correlations, whereas an irreducible hull presents M by generators 
and cogenerators. 

3. Computing with irreducible hulls 

Given a Q-graded module M in the usual way, via generators and relations, this 
section computes an irreducible hull M W as well as the cokernel of this inclusion. 

Calculating an effective irreducible hull of M is, by definition, equivalent to calcu- 
lating an irreducible decomposition of M. Thinking of the case M = k[Q]/I for a 
monomial ideal /, this procedure is polyhedral in nature: it writes the set of mono- 
mials outside of J as a union of convex polyhedral regions whose facets are parallel 
to those of Q. The algorithm for computing an effective irreducible hull M ^ W, 
culminating in Proposition 13. 7( does not require Q to be saturated. 

Computing the cokernel, however, is strictly easier for saturated semigroups. The 
main point is the computation of generators for irreducible ideals. For saturated 
semigroups this is Proposition 13. 141 The harder unsaturated case, in Proposition 13. 161 
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relies on the computation of irreducible ideals over its saturation. To highlight the 
simplification in the saturated case, we state the main result of Sections 13 . 2H3 .31 here . 

Proposition 3.1. Generators and relations for M and W/M are algorithmically com- 
putable from an effective irreducible hull M ^ W over any affine semigroup ring k[Q]. 

Proof. Generators for W are already given, and relations for W constitute a direct 
sum of irreducible ideals calculated as in Proposition 13.141 for saturated semigroups, 
and Proposition 13 . 161 in general. Since M is specified by its generators as a submodule 
of W, the current proposition reduces to calculating submodules and quotients of 
modules presented by generators and relations. □ 

3.1. Effective irreducible hulls from generators and relations. This subsection 
does not require the affine semigroup Q to be saturated. The next two results make 
Algorithm 13.61 possible to state and easier to read. The notation {yi, . . . ,yj) means 
'the k[(5]-submodule generated by the elements yi, ■ ■ . ,yj in their ambient module', 
and (0 -.M Pp) is the submodule of M annihilated by Pp. 

Lemma 3.2. Suppose F has minimal dimension among faces of Q such that Pp is 
associated to M. Then the natural map (0 :m Pp) to its localization (0 -.m Pf)\^F] 
along F is an inclusion. Furthermore, we can find algorithmically a set B G {0 :m -Pf) 
of homogeneous elements that consitute a k[ZF] -6aszs for (0 :m Pf)[^P]- 

Proof. The k[(5]-module (0 :m Pf) is naturally a torsion-free Ik[F]-module, by min- 
imality of dimF. Therefore (0 :m Pf) includes into its localization along F, which 
must be a free k[ZF]-module. Now use the following algorithm. □ 

Algorithm 3.3 (for Lemma (3. 2|) . Choose any element of (0 :m -Pf) as the first basis 
vector yi G B. Having chosen yj, let yj+i be any element of (0 -.m Pf) whose image in 
(0 :m Pf)/ {Ui^ ■ ■ ■ ,yj) generates a submodule of Krull dimension dimF (equivalently, 
the image of yj+i has annihilator Pp). The algorithm terminates when the Krull 
dimension of the quotient (0 :a/ Pp)/{yi, . . . ,yj) is strictly less than dimF. 

Lemma 3.4. In the situation of Lemma IJ.M the scalar factor on the (monomial) 
coefficient ofy&B in the unique k[7jF]-linear combination of elements in B equaling 
any fixed element 2 G (0 :m Pp) can be computed algorithmically. 

We present the proof as an algorithm. 

Algorithm 3.5 (for Lemma . Let B(z) = {y E B \ deg(y) = deg(z) (modZF)}. 
The coefficient of y in z is zero if y ^ P{z)- Otherwise, find elements a and {ay \ 
y G B{z)} in the face F such that a + deg(2;) = ay + deg{y) for all y G B{z). By 
construction, {x"^ ■ y \ y E B{z)} is a k-basis for the degree a + deg(z) piece of 
(0 -.M Pf), and standard methods allow us to calculate the syzygy with x"^ ■ z. 

Write TpN := Tp^N = (0 :iv Pp") for the set of elements in N annihilated by all 
high powers of Pp. 

Algorithm 3.6. 

INPUT Q-graded module M given by a generating set G C M and relations 
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OUTPUT effective irreducible hull W of M with effective vector set A indexed by G 



INITIALIZE 

W 



M 



({}, {}), the empty effective datum for the irreducible sum 
for all g E G; here () is the effective vector of length zero in W 

1 

DEFINE (Fi, . . . , Fg) := an ordering of the faces of Q with dim(Fj) < dim(Fj+i) 
WHILE i < s DO 

DEFINE F := Fi 

B := k[ZF]-basis for (0 -.n Pf)[^F], as in Algorithmic 
WHILE y E B and G G do 

IF (0 Pp) 7^ in some degree ayg = deg(y) (mod ZF) 



THEN Xyg := scalar coefficient of y on x 



g, as in 



ELSE A, 



Algorithm EH 





'yg ■ 

END IF-THEN-ELSE 
END WHILE-DO 

REDEFINE A^, := Concatenation of the two vectors \g and {Xyg)yeB, for g E G 
W := W ® {^B copies of F, Z^'-degrees of vectors in B) 

N := M/TpM 
i:= i + 1 
END WHILE-DO 

OUTPUT W along with A = {Agj^gG. where Xg is in degree deg{g) 

Proposition 3.7. Algorithm \3.b\ outputs an effective irreducible hull of M , using 
generators and relations for M as input. 

Proof. We must show that the homomorphism M ^ W determined by G and A is 
well-defined and injective. More precisely: monomial combinations z of the generators 
of M are zero if and only if the corresponding monomial combinations z\ of the A^ are 
zero in W] here, \g represents not a data structure but an element of W . 

The combination z is nonzero in M if and only if the submodule {z) C M generated 
by z has an associated prime. The associated prime is F := Fj if and only if the image 
of {z) in the succesive quotient A^ = M/Vp._^M intersects (0 itv Pp) nontrivially 
(this in particular imlplies that (0 -.^ Pp) is nonzero, so Pp is associated to M). 
This nontriviality of {z) fl (0 :iv Pp) is equivalent to having at least one of the terms 
monomial- (7 appearing in z be nonzero in the same (0 -.n Pp), because i? is a basis 
for (0 :n Pf)[^F]- Finally, monomial- (7 is nonzero precisely when the corresponding 
element monomial- Ag has nonzero coefficient in the appropriate summand of W. □ 

Remark 3.8. Some alterations to Algorithm 13.61 mav improve its running time. 
1. It is possible to avoid taking the successive quotients N/TpM at the redefine 
step. These quotients are designed to make Lemmas l3. 21 and IT^ applv. as well as 
to make A^ successively simpler. However, the cost of taking these quotients may 
not be worth it, since the final sentence of Lemma IT^holds even if F doesn't have 
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minimal dimension (so (0 -.m Pf) doesn't include into its localization along F). 
In fact, both of Algorithms 13.31 and 13.51 still work in this more general setting. 

2. Of the faces on the list (Fi, . . . , F^), only those associated to M need to be tested. 
If desired, these faces can be detected using homological methods. 

3. Instead of computing and working with (0 :m Pp) for each face separately, one 
could work with the modules (0 : m h) for each c, where Ic is the intersection of 
all primes Pp for faces F of dimension c. 

3.2. Generators and relations from irreducible hulls: saturated case. In this 
subsection we assume that Q is saturated. Our goal is to compute relations on the 
generators for M that come as part of an effective irreducible hull M W. As we 
shall see in the proof of Proposition 13.11 the computation essentially reduces to the 
case where M = W is an indecomposable irreducible sum W, so we are to determine 
the kernel of the surjection k[Q] —>■ W. More explicitly, given a face F and a degree 
a G Q, we must find generators of 

(1) W := k{Q\{a + F-Q)} 

as an ideal in k[Q]. 

Since Q is saturated, there is a unique minimal set of oriented hyperplanes inside Z'^ 
whose closed positive half-spaces in Z*^ have intersection equal to Q. The map sending 
H HnQ gives a bijection from these hyperplanes to the facets of Q. Denote by 
the closed positive half- space determined by an oriented hyperplane H, and by 
the open positive half-space. Thus is the complement of —H^ but can also be 
characterized as the lattice distance 1 translate of if+ in the positive direction. 

Lemma 3.9. Given any face F of Q and any element a & Q, 

Q\{a + F-Q) = [jia + Hl)nQ. 

Proof. We have a + F — Q = fl^Di^ ~ ^+ because Q is saturated (recall F — Q = 
-{Q + ZF)). Thus Z'^\{a + F -Q) = [j^^p a + H^. Now intersect with Q. □ 

Lemma [3.91 reduces the computation of generators for as in (0) to the case where 
F is itself a facet, at least when Q is saturated. The next algorithm and two lemmas 
cover this case by producing some rational polytopes whose integer points do the job. 
For notation, M+F denotes the real cone generated by F in M'^ = M (8) Z"', and 
denotes the real span of a hyperplane H . Also, by a Q-set we mean a subset of Z"' 
closed under addition by elements of Q. A set G of vectors in 7/ generates a Q-set 
T if T = G + Q. 

Lemma 3.10. Let Gq he the zonotope that is the Minkowski sum of all primitive 
integer vectors along rays of Q. Then, for all a G , the lattice points in a + Gq 
generate {a + R+Q) CiTj'^ as a Q-set. 

Proof. Let (3 he a. lattice point in a + M+Q. If there is no primitive integer vector p 
along a ray of Q such that (3 — p still lies in a + M+Q, then (3 E a + Gq. □ 

Algorithm 3.11. 
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F 

A 

INITIALIZE B 



INPUT Q := a saturated semigroup 

H := one of the hyperplanes bounding Q 

a e Q 

OUTPUT finite set 5 C Q such that the ideal (x^ | 6 G B) equals k{(a + n Q} 
DEFINE G := the polytope Gq in Lemma ITTUl 
H nQ, a facet of Q 

the set of faces of Q intersecting F only at G Q 
{}, the empty subset of Q 

WHILE D e A DO 

DEFINE Bf) := lattice points in Minkowski sum ((a + RH) n M+D) + G 
REDEFINE B := B U Bd 
NEXT D 
END WHILE-DO 
OUTPUT B 

Lemma 3.12. Alaorithm \3.11\ computes generators for the ideal Ik:{(a + H°^_) fl Q}. 

Proof. Suppose 6 G (a + if^) fl Q. The intersection {h + Mif) fl M+Q is a polyhedron 
whose bounded faces are precisely the polytopes (h + M.H) fl M.^D for D E A, and 
whose recession cone is M+F. Therefore h E h' + M+F for some real vector h' G 
(6 + Mil) n R+D and some face D G A. Moreover, h' lies in 6" + M_|_D for some real 
vector h" G (a + RiJ) fl M+D. Consequently, h lies in 6" + R+{D + F), and therefore in 
6" + M+Q. Now x'' lies in the k[(5]-module generated by klS^}, by definition of G. □ 

Remark 3.13. Some alterations to Algorithm I3.11l mav improve its running time. 

1. Instead of computing just one polytope G = Gq and Minkowski summing it to 
define every Bd, we could define Bd with Gd+f in place of G, for each face 
Z) G A. This might reduce the number of lattice points in B dramatically, but 
would require more computations as in Lemma f3.1Ul 

2. Restricting to the maximal elements in A will speed things up. 

Let us summarize the above algorithm and three lemmas. (See Section [7| for is- 
sues concerning the output of the algorithm in the following proposition, and post- 
processing for the purpose of reducing its complexity.) 

Proposition 3.14. Generators of the irreducible ideal W = ker(k[(5] —>■ W) are algo- 
rithmically computable using as input an indecomposable effective irreducible sum W 
over a normal semigroup ring k[Q]. 

Proof. Apply Algorithm 13. 1 II to each of the sets (a -|- H^) n Q in Lemma [3.91 □ 

3.3. Generators and relations from irreducible hulls: unsaturated case. Now 

we return to the general case, where Q need not be saturated, and denote by Q^^^ 
the saturation of Q. The basic idea for computing generators of irreducible ideals 
in k[Q] is to intersect (as k[(5]-modules) the submodule k[Q] C k[(5'^^*] with the ideal 
W C k[(5*^'^*] output in the saturated case. Proposition 13.141 Then it remains to find 
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the appropriate F-primary component of as a k[(5]-niodule, where F is the unique 
face of dimension dim(iy) associated to W (as a k;[(5]-module). 

Every module in Algorithm 13. 151 is to be considered as a k[(5]-module — even those 
generated as k[(5'^'^*]-modules. Thus F is always a face of Q, and we consider F — Q as 
opposed to F — Q^^^. Note, however, that Ik{F — Q^^^} does equal the corresponding 
injective over Ik:[(5''^*], even though F is a face of Q; subtracting Q^^^ automatically 
saturates F. 



Algorithm 3.15. 



INPUT 


Q 






F 






a 


G 


OUTPUT 


B 


C 


DEFINE 


V 






V 






w 






I 




INITIALIZE 


B 






w 




WHILE 


(0: 


w 


DEFINE 



a semigroup, not necessarily saturated 
a face of Q 

Q 

Q such that (x^ | 6 G B) equals the ideal k{Q \ (a + F - Q)} in k[Q] 

kjQ^'^* \ (a + F - Q""""^)}, an indecomposable irreducible over kfQ'*''*] 

the kernel of k[Q'''*] V output by Proposition ITT41 

V^nk[(5], the intersection taken inside kfQ'*'^''] 

f] {Pd I -D is a facet of F}, an ideal in k[Q] 

degrees of the elements generating W 

k[Q]/W 

Pp) has a generator in some degree ^ a (mod ZF) DO 



REDEFINE 

DEFINE 
REDEFINE 



END 
OUTPUT 



G 
B_ 

W 
G' 
B 

w 

WHILE-DO 
B 



generators for (0 Pp) that lie in degrees ^ a (mod ZF) 
B U degrees of the elements in G 

W/G 

generators for TiW 

B U degrees of the elements in G' 

W/G' 



Proposition 3.16. Alaorithm VJ.l^ outputs generating degrees for k{Q\{a + F — Q)} . 

Proof. The module W gets initialized as a quotient of k[Q'^'^*] with dimension dim(F) 
as a k[(5]-module. This much holds by the saturated version Proposition 13. 141 applied 
to V, and the preservation of dimension |Eis95| Proposition 9.2] for the module-finite 
ring extension k[Q] C k[(5'''^*] applied to V. One part of the output is clear: the set 
(x* I 6 G -B) generates the kernel of the map k[Q] W at every stage in the algorithm. 
The question is whether W is the claimed indecomposable irreducible sum. 

In the first REDEFINE step, the annihilator of x" G remains Pp. Indeed, any 
element killed by Pp that generates a submodule containing a nonzero element in 
degree a must itself have degree congruent to a (modZF). The second redefine step 
only kills elements with annihilators strictly larger than that of x*^; such elements can't 
generate submodules containing x". Therefore, W has only one associated prime Pp 
after each loop of while-do, by dimension considerations. 
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When the loop terminates, the locahzation (0 :pp Pp)[ZF] along F is indecom- 
posable, being isomorphic to k{a + ZF}. It follows that the kernel of the surjec- 
tion k.[Q] — i> is an irreducible ideal |Vas98t Proposition 3.1.7]. We are done by 
Lemma \2A\ because k{a + F — Q}q is the only indecomposable irreducible sum for 
which the annihilator of x'^ is Pp- □ 

Remark 3.17. Some alterations to Algorithm 13 . 1 51 mav improve its efficiency. 

1. The step W := W /TiW need not occur until the very last step before output. 
Its current placement is designed to speed the computation by simplifying W in 
each loop, but the cost of taking the colon may not make up for it. Instead, the 
end of the algorithm can be replaced by: 

WHILE (0 :p57 Pp) has rank strictly larger than 1 over k[F] DO 

DEFINE G := generators for (0 -.jy Pp) lying in degrees ^ a (modZF) 
REDEFINE B := B U degrees of the elements in G 

W := W/G 

END WHILE-DO 

REDEFINE B := B U degrees of the generators of TjW 

OUTPUT B 

2. As in Remark |3.8^ it is not necessary to compute all of (0 ipp Pp) in the WHILE- 
DO loop. It suffices instead to let G be a basis for (0 :pp Pp)[ZF]. This remark 
also holds for the reworked while-do loop in the previous item. 

3. The set B can become rather redundant. Since the machine will have to keep a 
presentation of W in memory, the algorithm could simply spit out the relations 
defining ly as a k[Q]-module at the very end, without keeping track of B at all. 

4. Computing injective resolutions 

In this section the semigroup Q is not required to be saturated. Our goal is the 
main result (Theorem 14. 7|) in the first half of the paper: an algorithm to compute 
injective resolutions of finitely generated modules over k[Q], in the Z'^-graded setting. 
That is, given generators and relations for a finitely generated Z'^-graded module M, 
we will compute an exact sequence ^ M ^ J° ^ ^ ■ ■ ■ in which J* is a Z"^- 
graded injective module for each i. Of course, we shall only say how to calculate up to 
some specified cohomological degree, as injective resolutions usually do not terminate. 
This will not pose a problem for our subsequent computation in Section [7| of local 
cohomology, which vanishes past cohomological degree d + 1 anyway. 

The upshot is to reduce the computation of injective resolutions to finding irreducible 
hulls of finitely generated Q-graded modules and computing their cokernels, which we 
have already done in Sectional 

The data structures we employ for Z^'-graded injective resolutions are the matrices 
we introduce in the next definition. 

Definition 4.1. A monomial matrix is a matrix of constants Xqp along with 

1. a vector aq G Z'^ and a face Fg & Q for each row, and 

2. a vector Op G Z'^ and a face Fp & Q for each column 
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such that Xqp = unless Fp C Fq and Op G aq + Fq — Q . 

These monomial matrices generalize those in |MilOOj . which were for Q = N'^. 

To any monomial matrix we can associate a map J i— > J' of injective modules in the 
following manner. Each row and column label gives the data of an indecomposable 
injective; we think of the row labels as giving summands of J and the column labels 
as giving summands of J'. To give a map from J to J' is thus the same as giving a 
matrix of maps from the row indecomposables to the column indecomposables. Such 
a map k{aq + Fq — Q} i— > k{ap + Fp — Q} is necessarily zero unless Fp C Fq and 
ap & aq + Fq — Q. In the latter case it is determined by a single scalar Xqp. Hence 

••• Fp ... 



Oip 




is a monomial matrix representing a map 

0k{a, + F,-g} ^ 0k{ap + Fp-Q}. 

q p 

The component k{aq + Fq ~ Q} i-^ k{ap + Fp — Q} of this homomorphism takes 
to Agpx" for aA\ a & ap + Fp — Q , and is zero elsewhere. 

Note that in degree a, the map Ja ^ J'a given by a monomial matrix is obtained by 
deleting the rows and columns labeled by ap, Fp such that a does not lie in ap + Fp — Q. 
(This corresponds to ignoring those summands of J and J' not supported at a.) 
Ignoring the labels on what remains gives us a matrix with entries in k, which defines 
the k- vector space map Jq ^ J'^- 

Two monomial matrices represent the same map of injectives (with given decom- 
positions into direct sums of indecomposable injectives) if and only if (i) their scalar 
entries are equal, (ii) the corresponding faces Fr are equal, where r = p,q, and (iii) the 
corresponding vectors Or are congruent modulo ZF^. 

Rather than compute directly with cumbersome, infinitely generated injectives, it 
is more convenient to approximate injective resolutions using irreducible sums. 

Definition 4.2. An irreducible resolution of a Q-graded module M is an exact 
sequence — M — — > — > ■ ■ ■ in which each is an irreducible sum. 

Irreducible resolutions are approximations to injective resolutions; indeed, the Q- 
graded part of any injective resolution is an irreducible resolution jMil02t Theorem 2.4]. 
In particular, monomial matrices just as well represent homomorphisms of irreducible 
sums, as long as the degree labels Uq and ap all can be chosen to lie in Q. The (appar- 
ent) advantage to irreducible resolutions over injective resolutions is their finiteness. 

Corollary 4.3. For any finitely generated Q - graded k[Q]-module M, Propositions \S. 1\ 
and \3. 7| inductively compute a minimal irreducible resolution W of M algorithmically. 
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Proof. Minimal irreducible resolutions have finite length (that is, they vanish in all 
sufficiently high cohomological degrees) by [Mil02l Theorem 2.4]. The computability 
therefore follows from Propositions 13.11 and 13.71 bv induction on the highest cohomo- 
logical degree required. □ 

The next result demonstrates the precise manner in which irreducible resolutions 
approximate injective resolutions for computational purposes. 

Proposition 4.4. Let M be a finitely generated module with minimal injective reso- 
lution J' and minimal irreducible resolution W . Suppose that every indecomposable 
summand in the first n cohomological degrees of J' has nonzero Q-graded part. Then 
M is Q-graded, and the data contained in the first n stages of W constitute a finite 
data structure for the first n cohomological degrees of J' . 

Proof. Every map in J' can be expressed using the finite data of a monomial matrix, 
and this data can be read immediately off the maps in ly* . □ 

If we can algorithmically determine a Z'^-graded shift of M so that the hypotheses 
of Proposition 14.41 are satisfied, then we can compute the minimal injective resolution 
of M up to cohomological degree n. This task requires a lemma, in which m denotes 
the maximal ideal P{o} generated by nonunit monomials in k[(5]. 

Lemma 4.5. Let J' be a minimal injective resolution of a finitely generated module M , 
and F a face of Q. If every indecomposable summand of ji+'^-dim(F) nonzero 
Q-graded part, then every indecomposable summand of isomorphic to a U^-graded 
shift of lk{F — Q} has nonzero Q-graded part. 

Proof. |HM03t Proposition 3.5], in the special case of an affine semigroup ring. □ 

Every indecomposable summand of T^J^ is a shift Ikja — Q} of k;{— Q}- Such 
an indecomposable injective has nonzero Q-graded part if and only if a G Our 
final lemma in this section describes the (standard) way to calculate the shifts a 
appearing in Pm^-'. The number /i-^'"(M) of shifts kjo; — Q} appearing as summands 
in cohomological degree j of the minimal injective resolution of M is called the j*^ 
Bass number of M in degree a. 

Lemma 4.6. Let T. be a free resolution of the residue field k. The Bass number 
/i-''"(M) is effectively computable as the k-vector space dimension of {T.., M)^. 

Proof. This expression of Bass numbers as dimensions (over k) of Ext modules is 
standard; see |BH93l Chapter 3]. The computability follows because we can calculate 
free resolutions, homomorphisms, and homology over k[Q]. □ 

Now we come to our central result. For notation, M(— a) denotes the Z'^-graded 
shift of M up by a, so that M{—a)b = Mb^a- 

Theorem 4.7. Fix a finitely generated k[Q]-module M and an integer i. There is 
an algorithmically computable a & Q for which Propositions IJ.il and \3. 7| inductively 
compute the minimal injective resolution of M{—a) through cohomological degree i + 
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Proof. After using Lemma to compute the Bass numbers of M up to cohomological 
degree i + 1 + dim(M), choose a so that the corresponding Bass numbers of M(— a) 
have Z'^-graded degrees lying in Q. At this point, M(— a) satisfies the hypotheses of 
Proposition 14.41 with ?7, = i + 1, by Lemma [4.51 Now apply Corollary 14.31 □ 

5. Sector partitions from injectives 

We turn now to sector partitions, for which we assume henceforth that the affine 
semigroup Q is saturated. As a prerequisite to producing sector partitions of local 
cohomology modules, we demonstrate in this section that injective modules admit 
sector partitions, as does the homology of any complex of injective modules. 

Proposition 5.1. Suppose J = Ji is an injective module decomposed into sum- 
mands Ji = k{ai + Fi — Q}. For each subset A <Z {1, . . . ,r} define Sa to be the set 

= {aeZ'^l (Ji)^ ^kforieA} 

of all degrees in Z'^ such that the summands of J nonzero in that degree are precisely 
those indexed by A. The sets Sa canonically determine a sector partition S{J) h J. 

Proof. For each a G Z"', either {Ji)a = {0} or {Ji)a = k-x". Therefore S{J) is indeed a 
partition of Z'^. Now we must show that Sa is a finite union of polyhedra as in part 1 of 
Definition II. 21 The set a + F — Q oi degrees is the set of lattice points in a polyhedron 
of the desired form because the half-spaces whose intersection is a + F — Q are bounded 
by hyperplanes parallel to facets of Q, by definition. These hyperplanes divide Z'^ into 
finitely many disjoint regions (place the lattice points lying on each hyperplane in the 
region on the positive side of that hyperplane), each of which consists of the lattice 
points in a polyhedron of the desired form. Thus the complement 'Z'^\{a + F — Q) is 
the required kind of finite union. We conclude that Sa is a finite union of regions, each 
of which is an intersection of r polyhedral regions — one from each of the summands J,. 

For each index set A such that Sa is nonempty, define C k*" to be the subspace 
spanned by the basis vectors Cj such that i & A. Then for each degree a in 5a, the 
map Ja Jsa required by part 2 of Definition II . 21 can be taken to equal the zero map 
on {Ji)a for i not in A, and the map sending x" to Cj on {Ji)a for i in A. 

To define the maps x'^^"'^-* for index sets A and B such that Sb — Sa is nonempty, 
as in part 3 of Definition 11.21 it suffices to define the image of for each i in A. We 
take x'^^~'^'*(ej) = Ci if i is in B, and x'^'^~'^-^(ej) = otherwise. Commutativity of the 
required diagram follows from the definition of the module structure on kjaj + Fj — Q}. 
Specifically, for a E Sa and (3 & Sb with P — a E Q, multiplication by x^~" takes x" 
to x^ in Ji for i ^ B, and takes x° to zero in Jj for i outside B. □ 

The sector partition in Proposition l5. ll descends to the cohomology H of any complex 
of injectives, via monomial matrices. The forthcoming sector partition of H is really 
determined canonically by J* (without its direct sum decomposition), even though 
the way we present things here makes it look like bases must be chosen. We chose this 
route because bases are good for computation, while uniqueness is immaterial. 
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Theorem 5.2. If H is a module that can be expressed as the (middle) homology of 
a complex J' : J' ^ J ^ J" in which all three modules are injective, or all three 
modules are flat, then there is a sector partition S{J') h H determined by J' . 

Proof. Choose direct sum decompositions to write 

r' r r" 

J' = J'i-, J = Ji-, J" = . 

i=l i=l 1=1 

Let $ and \1/ be the monomial matrices representing the maps J' ^ J and J J", 
respectively. The sectors in the sector partition iS(J' © J © J") h J' ® J ® J" are 
indexed by triples (A', A, A") of subsets of {1, . . . , r'}, {1, . . . , r}, {1, . . . , r"}, respec- 
tively, and automatically satisfy the polyhedrality condition in part 1 of Definition II .21 
by Proposition 15.11 We take S{J') to partition Z*^ into these sectors. 

For each triple {A', A, A") we have maps : Js^, Jsj^ and : Js^ — > J's^,, 
whose monomial matrices are defined by deleting: row i' of $ for i' not in A'; column i 
of $ and row z of \1/ for i not in A; and column i" of \1/ for i" not in A". Let 

(2) Hs,,,^,„ = ker«)/im($:^,)- 
For any a in Sa',a,A", we have a commutative diagram 

7' > 7 ^ J" 

(3) i $A i ^a" J- 



T/ -4' 7 ^ , 111 

. — ^ Jsa — ' -'s 



A" 



that induces the required isomorphism = Hs^, ^ ^„ . It is routine to check that the 
maps Hsj^i ^ ^„ — > -f^s^, ^ g„ induced from the corresponding maps on 

J A' ^ J and J A" 

commute with this isomorphism. □ 

Once we have Theorem 15. 2[ the only step remaining to prove Theorem 11.31 is to 
exhibit H]{M) as the homology of a complex of injectives. 

Remark 5.3. The results in this section hold just as well for fiat objects of A^, 
which are Matlis dual to injective objects and hence isomorphic to finite direct sums 
of modules of the form k{a + F + Q} for some a in T/' and some face F of Q \M.ii03\ 
Chapter 11]. For the proofs, simply apply Matlis duality to the results for injectives. 

6. Computing sector partitions 

Again letting Q be a saturated affine semigroup, the next task is actually computing 
the finitely many polyhedra whose lattice points comprise the sectors in the sector 
partition S{J) h J of an injective module. That is, we need to make Proposition 15.11 
and its proof into an algorithm. 

Since Q is saturated, there are unique primitive integer linear functionals ri, . . . , r„ 
taking Z*^ Z, one for each facet of Q, such that Q = HiLii'^* — 0} 
lattice points in the intersection of their positive half-spaces. The degrees on which 
indecomposable injectives are supported can be expressed in terms of these linear 
functionals, via the following identity: 

(4) a + F~Q = {f3 e Z'^ \ Ti{l3) < Ti{a) whenever Ti{F) = 0}. 
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In other words, F — Q is the intersection of the negative half-spaces for those function- 
al Tj vanishing on F, and a + F — Q is simply a translate. By convention, we use the 
notation rj(/9) < oo to mean that there is no restriction on the value of Ti{f3). This 
allows a notation Tp{a) G (Z U oo)"' for the vector whose i^^ coordinate satisfies 

1 oo otherwise. 

The point is that a vector /5 G Z'^ lies in a + F — Q if and only if t{i3) < TF{a), where 

r(/3) = (n (/?),..., r„(/5)) 

and the '<' symbol denotes componentwise comparison. We shall use the correspond- 
ing definitions of Tp{a) and t{/3) for vectors G = M®Z°', so Tp{a) G (MUoo)'^. 
For the rest of this section, let 

r 

(5) J = with = k{aj + Fj-Q}, 

be an inject ive module, and define 

:= Tp^iaj) for j = 1, . . . ,r. 

Thus for i = 1, . . . ,n the vector has i^^ coordinate t- = r^?, {aj)i, which equals either 
Ti{aj) or oo, depending on whether Tj vanishes on Fj or not. Even without calculating 
the set S{J) algorithmically, the vectors specify the map from Z'^ to 5 (J), by 
definition. We record a precise version of this statement in the next lemma. 

Lemma 6.1. A degree a E Tj'^ lies in Sa if and only if A = |j G {1, . . . , r} 

It remains to ascertain which sets £14 of lattice points are nonempty, and to determine 
the pairs A,B for which we must compute a map x"^"^ : — > Jb- (The maps 
themselves, which are canonical, are constructed in the proof of Proposition 15.11 ) For 
each functional there is a permutation Wi of {1, ... , r} satisfying t^^^^ < ■ ■ ■ < t^'^'^\ 
To simplify notation, we write f/ instead of t^''^^\ Also, set f° = —00 and f[~^^ = 00. 

For fixed i, the parallel affine hyperplanes {rj = divide Z'^ into strips 

for i = 0, . . . ,r. At most r + 1 of these strips are nonempty, because some of the 
hyperplanes may coincide. Also, the last few of the f/ will equal 00; we interpret any 
strip where r/ = r^^^ = 00 as empty, and ignore it. 

Proposition 6.2. Let J be as in (0). For any fixed £1, G {0, ... , r}, the lattice 

points in the polyhedron 

n 

A(£i,...,4) := []{f3eR''\f!^+l<T,{P)<f^^^'} 

i=l 

all lie inside a single sector in S{J). The partition of Z'^ by the polyhedra A(£i, . . . , in) 
refines the partition of Z^ by the sectors in S{J). 

Proof. This follows from the definitions and (HJ, which uses that Q is saturated. □ 
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Proposition 16.21 makes way for an algorithm to compute the set of sectors. 
Algorithm 6.3. 

J = 0^=1 , an injective module over k[Q], with .P = k{aj + Fj — Q} 
the set iS( J) of sectors, each expressed as a list of polyhedra that partition it 
subsets of {1, . . . , r}, as in Lemma 16.11 
{}, the empty collection of subsets of {1, . . . , r} 
,£„ G {0,...,r} DO 
A(£i,...,C)7^0 

THEN DEFINE A := 0(A(£i, ... ,£„)) 
ELSE NEXT (4, . . . , £„) 



INPUT 
OUTPUT 
DEFINE 
INITIALIZE 
WHILE 



: 

A: = 
IF 



END 
OUTPUT 



END IF-THEN-ELSE 
IF Ae A 

THEN REDEFINE 
ELSE INITIALIZE 
REDEFINE 
END IF-THEN-ELSE 
NEXT 
WHILE-DO 

{SaIAg A} 



SA:=SAU{A{h,...,Q} 

SA:={A{i,,...,Q} 
A:=AU{A} 



by definition, and can easily be determined 



Note that is constant on A(£i, . . . , 
directly from the data (£i, . . . , 

Next comes the determination of which maps :k^^^ need computing. In the com- 
ing algorithm, we write A(£i, ...,£„)< if ...,£„) < . . . , 4) as 
vectors in (Z U oo)". Such notation is justified because A(£i, . . . ,in) ^ A.{£[, . . . , 
automatically implies that . . . ,£'^) — A(£i, ...,£„) fails to intersect Q. 

Algorithm 6.4. 

INPUT sectors Sa and Sb in S{J) from the output of Algorithm 16.31 
OUTPUT the truth value of: "there exist a e Sa and P e Sb with (3 — a e Q" 
INITIALIZE val := FALSE 

WHILE (Aa, Ab) e Sax Sb and val = 
IF AD B and Aa < Ab 

THEN DEFINE A^ - A^ 
ELSE NEXT {Aa, Ab) 
IF-THEN-ELSE 

Q n (Ab - Aa) ^ 

THEN REDEFINE val : = 
ELSE NEXT {Aa,Ab) 
IF-THEN-ELSE 



FALSE, DO 



the Minkowski sum of A^ and — A^ 



END 
IF 



TRUE 



END 
OUTPUT 



END 
WHILE-DO 

val 



The proof of correctness for Algorithm 16.41 is straightforward from the definitions, 
except for the first if-then-else procedure, which relies on Lemma f6. 51 below. Note 
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the non-necessity in Algorithm l6 . 4l of actually finding a witness in A^— for Sa ^ Sb] 
as we have seen in (j2)) and from the proof of Theorem 15.21 the natural map on 
cohomology is induced by taking submatrices of the monomial matrix, regardless of 
where the witnesses lie. 

Lemma 6.5. IfS{J) is as in Proposition \5.1l then QCi {Sb — Sa) ^ implies AD B. 

Proof, li a E Q and {Ji)a = 0, then {Ji)a+a = 0, so the set of summands nonzero in 
degree a + a can only be smaller. □ 

Unfortunately, Algorithm 16.41 is necessary, because A^ — A^ 7^ need not always 
hold when A^ < A^, as the example to come shortly demonstrates. It does seem, 
however, that the offending pairs of polytopes are usually "small". For instance, we 
know of no examples where the lattice points in either polytope affinely span Z'^. 

Example 6.6. Let Q be the subsemigroup of generated by (2, 0), (1, 1), and (0, 2). 
Name the faces of Q as 0, X, Y, Q, and set Ep = F — Q. Let 

J = k{(0,0) + Eo} © k{(0,l) + Ex} © k{(0,0) + Ey} 

© k{(0, -1) + Ex} © k{(-2, 0) + Ey}, 

with the summands labeled in order as Ji, . . . , J5. Letting X be facet number 1 and 

Y be facet number 2, the arrays t- and f/ look like 

rrl\ /O 1 00 -1 00 X ^ /-oo -l O l oo oo ooX 

\T^J V0oo0oo-2y \f^J V-oo-2 00ooooooy 

The sectors 5'{i,2,3} and S{2,3} contain one polytope each, and both of these polytopes 

contain exactly one lattice point. Specifically, identifying the sector, the polytope, and 

the lattice point, we have 

^{1,2,3} = A(-l, -2) = (0, 0) and 8^2,3} = A(0, -2) = (-1, 1). 

Now A(— 1, —2) < A(0, —2), but subtracting the vector in 5'{i,2,3} from the one in S{2,3} 
yields (—1, 1), which does not lie in the semigroup Q. 

Remark 6.7. The notion of sector partition ought to have a refinement that takes 
into account the various kinds of failures of saturation for arbitrary affine semigroup. 
The resulting notion would produce sector partitions for the cohomology of complexes 
of injectives over nonnormal affine semigroup rings. The failures of saturation fall into 
two categories: the geometric kind, arising from polyhedral "holes" in the semigroup 
(as compared with its saturation), and the arithmetic kind, arising from finite-index 
sublattices generated by faces. Even in the case where arithmetic failure is absent, 
however, we do not know how to bound the sizes and shapes of the "holes" sufficiently 
to carry out an analysis such as the one producing the algorithms above. 

7. Computing local cohomology with monomial support 

Still assuming that Q is sturated, we have now finally developed enough tools to 
prove the main theorem on local cohomology with monomial support, namely Theo- 
rem from the Introduction. 
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Proof of Theorem II .31 Take i = d in Theorem Wli and let J' {-a) be the min- 
imal Z'^-graded injective resolution computed there. Then J' is an algorithmically 
computed injective resolution of M. By definition, H]{M) is the middle cohomology 
of the complex F/J*"^ — > F/J* F/J*"*"^, where F/J-' is the direct sum of all inde- 
composable summands of whose unique associated prime contains /. Having now 
expressed H]{M) as the cohomology of an effectively computed complex of injectives, 
Theorem 15.21 savs that H}{M) has a sector partition. The set of sectors in part 1 of 
Definition 11.21 is computed by Algorithm 16.31 The vector spaces in part 2 of Defini- 
tion 11.21 are specified in Q from the proof of Theorem 15.21 and naturally determine 
the maps in part 3 of Definition 11.21 given the computation in Algorithm 16.41 D 

Now we turn to issues of complexity. There is little sense in completing a formal 
complexity analysis of all of the algorithms presented in this paper, as they involve 
Grobner basis computation, which is doubly-exponential from a worst-case perspec- 
tive. However, it is worth mentioning where the complexity in our algorithms comes 
from, up to a factor arising from the complexity of Grobner basis computation, since 
Grobner basis computations are often more efficient than expected. The purpose of 
what follows, therefore, is to assure the reader that our algorithms have not amplified 
the faux-doubly-exponential complexity of Grobner bases with some "honest" expo- 
nential complexity. 

Let us assume that the dimension d is fixed, and analyze the complexity of computing 
all of the local cohomology of a finitely generated module M supported on a fixed 
monomial ideal / over a normal semigroup ring k[Q]. This computation involves 
all of the algorithms in the paper except the one in Section 13.31 (The complexity 
of Algorithm 13.151 above and beyond Algorithm 13.61 is only about as bad as that of 
kfQ'^^^J/kfQ] as a k[(5]-module, anyway.) 

In Algorithm 13.61 the only non-Grobner contribution to the running time comes 
from the number of basis elements constructed (see Remark 13.81 2. which can be used 
to ensure that we only check faces of Q giving rise to basis elements). This number is by 
definition a Bass number of M. Thus, up to Grobner basis computation. Algorithm 13. 61 
is only as complex as its output. 

Next we consider the algorithm in Proposition 13. 14l The algorithm works by taking 
the union (over a set of facets of Q) of ideals output by Algorithm 13.111 The output 
presents the generators of each such ideal as the lattice points in a union of polytopes 
having the form ((a-|-Mif)nR+D) +G, where D is a face of Q. The computation of each 
such polytope is by standard techniques to intersect polyhedra and take Minkowski 
sums with the fixed zonotope G. Hence, up to factors coming from the number of 
facets of Q and from standard procedures, we need only bound 

1. the number of polytopes output by Algorithm 13.111 and 

2. the number of lattice points in each such polytope. 

The former is polynomial in the number of facets of Q by Remark 13.131 2. The latter 
is polynomial in the input vector a G Q by the piecewise polynomiality of the lattice 
point enumeration function of (a + Mif) flM+D as a function of a (McM77j . along with 
the fact that G is fixed. Actually computing the set of lattice points in each polytope 
can be accomplished using the efficient algorithms of Barvinok and Woods |BW02j . 
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Remark 7.1. We need to do Grobner basis computations with the irreducible ideals W 
output by the algorithm in Proposition 13.141 This means that, for our purposes, the 
short rational generating functions output by the algorithms of |BW02j do not suffice: 
we actually require the list of lattice point explicitly, to get a generating set of IV as a 
list of monomials. Thus the short generating functions must be expanded. To reduce 
complexity, the short generating functions can be post-processed using the methods of 
|BW02j to yield short generating functions for the minimal generators of the ideals W 
in question. Then we can expand only these "minimal" short generating functions. 

The remaining contributions to the complexity of our local cohomology computation 
come from Algorithm 16.31 which computes the sets of polytopes whose disjoint unions 
constitute the sectors, and Algorithm 16.41 The latter is quadratic in the output of 
Algorithm 16.31 times a factor coming from the Minkowski sum operations and the 
decision procedure for whether each such sum contains a lattice point after intersecting 
with Q. Therefore it remains only to analyze Algorithm 16.31 

Proposition 7.2. The number of polyhedra arising in Alaorithm \6.^ is polynomial in 
the Bass numbers of M and the number of facets of Q. 

Proof. Each Bass number of M represents an indecomposable injective module whose 
bounding hyperplanes subdivide into a number of regions. Consider the subdivision 
of obtained by taking simultaneously all of the hyperplanes corresponding all of 
the Bass numbers of M. The number of hyperplanes contributed by each Bass number 
is at most the number of facets of Q, so the total number of hyperplanes is at most 
the number of facets of Q times the sum of the contributing Bass numbers. It is 
well known (and follows by induction on n and the dimension d) that n hyperplanes 
subdivide into a number of regions that is a polynomial in n of degree d. □ 

This proof shows that the number of polyhedra is exponential in the dimension. 
Exponential growth as a function of dimension also occurs in the analysis before Re- 
mark [THl where we apply |McM77j . 

Remark 7.3. A large number of rational polyhedra arise in the course of computing 
local cohomology modules. When the identification of all the lattice points in these 
polyhedra is necessary, the complexity of this task should be drastically reduced by 
the fact that most of these polyhedra have facets parallel to those of Q itself. Results 
such as those in |BV97j could be helpful along these lines. 

Acknowledgments. We wish to thank Gennady Lyubeznik, who motivated us to 
make our previously abstract methods algorithmic. Both authors were partially sup- 
ported by the National Science Foundation. 

References 

[BH93] Winfried Bruns and Jiirgen Herzog, Cohen-Macaulay rings, Cambridge Studies in Ad- 
vanced Mathematics, voL 39, Cambridge University Press, Cambridge, 1993. 

[BV97] Michel Brion and Michele Vergne, Residue formulae, vector partition functions and lattice 
points in rational polytopes, J. Amer. Math. Soc. 10 (1997), no. 4, 797-833. 



22 



DAVID HELM AND EZRA MILLER 



[BW02] Alexander Barvinok and Kevin Woods, Short rational generating functions for lattice point 

problems, J. Amcr. Math. Soc, to appear. arXiv:math.CO/0211146 
[Eis95] David Eisenbud, Commutative algebra, with a view toward algebraic geometry, Graduate 

Texts in Mathematics, vol. 150, Springer- Verlag, New York, 1995, first printing. 
[EMSOO] David Eisenbud, Mircca Mustata, and Mike Stillman, Cohomology on toric varieties and 

local cohomology with monomial supports, J. Symbolic Comput. 29 (2000), no. 4-5, 583-600. 
[Har70] Robin Hartshorne, Affine duality and cofiniteness, Invent. Math. 9 (1969/1970), 145-164. 
[HMOS] David Helm and Ezra Miller, Bass numbers of semigroup- graded local cohomology, Pacific 

J. Math. 209 (2003), no. 1, 41-66. 
[HS93] Craig Huneke and Rodney Sharp, Bass numbers of local cohomology modules. Trans. Amer. 

Math Soc. 339 (1993), 765-779. 
[Lyu93] Gennady Lyubeznik, Finiteness properties of local cohomology modules ( an application of 

D-modules to commutative algebra). Invent. Math. 113 (1993), no. 1, 41-55. 
[McM77] P. McMuUen, Valuations and Euler-type relations on certain classes of convex polytopes, 

Proc. London Math. Soc. (3) 35 (1977), no. 1, 113-135. 
[Mil98] Ezra Miller, Alexander duality for monomial ideals and their resolutions, math.AG/9812095, 

1998. 

[MilOO] Ezra Miller, The Alexander duality functors and local duality with monomial support, J. 
Algebra 231 (2000), 180-234. 

[Mil02] Ezra Miller, Cohen- Macaulay quotients of normal semigroup rings via irreducible resolu- 
tions. Math. Res. Lett. 9 (2002), no. 1, 117-128. 

[MS03] Ezra Miller and Bernd Sturmfels, Combinatorial commutative algebra, version of 30 July 
2003 (currently available at http://math.umn.edu/~ezra), in progress. 

[MusOO] Mircea Mustafa, Local cohomology at monomial ideals, J. Symbolic Comput. 29 (2000), 
no. 4-5, 709-720. 

[Ter99] Naoki Terai, Local cohomology modules with respect to monomial ideals, preprint, 1999. 

[Vas98] Wolmer V. Vasconcelos, Computational methods inn commutative algebra and algebraic 
geometry, Algorithms and Computation in Mathematics, vol. 2, Springer- Verlag, Berlin, 
1998, With chapters by David Eisenbud, Daniel R. Grayson, Jiirgen Herzog and Michael 
Stillman. 

[Wal99] Uli Walther, Algorithmic computation of local cohomology modules and the cohomological 

dimension of algebraic varieties, J. Pure Appl. Algebra 139 (1999), 303-321. 
[YanOl] Kohji Yanagawa, Sheaves on finite posets and modules over normal semigroup rings, J. 

Pure Appl. Algebra 161 (2001), no. 3, 341-366. 
[Yan02] Kohji Yanagawa, Squarefree m,odules and, local cohom,ology modules at monomial ideals, 

Local cohomology and its apphcations (Guanajuato, 1999), Lecture Notes in Pure and 

Appl. Math., vol. 226, Dekker, New York, 2002, pp. 207-231. 

Mathematics Department, UC Berkeley, Berkeley, CA 
E-mail address: dheliii@iiiath.harvard.edu 

Mathematical Sciences Research Institute, Berkeley, CA 
E-mail address: ezra@math.umn.edu 



